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•^ , Abstract 

,^ ■ An exactly solvable model is used to investigate the assumptions behind color 



transparency. 



1 INTRODUCTION 

Color transparency (CT) is the anomalously high transparency of the nucleus to nucle- 
ons in quasi-elastic high-momentum transfer nuclear processes, measured with resolution 
good enough to insure that no "extra" pions are produced. This means that the ab- 
sorptive nuclear optical potential representing initial and final state interactions plays no 
role in such reactions. Color transparency is under active investigation with experiments 
performed at BNL |jl| , SLAC 0| and ongoing work at BNL and experiments proposed 
at CEBAF. Further references can be found in the review Ref. HI. 



The existence of color transparency depends upon three assumptions [^, |[ : 

(i) A small wave packet is formed in a high momentum transfer reaction. This wave 
packet is sometimes dubbed a point like configuration (PLC). 

(ii) The interactions between a small color neutral wave packet and the nucleus are 
suppressed. 

(iii) The wave packet escapes the nucleus before expanding. The expansion time r is 
typically stated as r ~ tqP/M where tq is a time in the rest frame (expected to be about 
1 fm) needed for a small system to expand and P/M is a time dilation factor. For large 
enough momentum, P, r is large. 

Each of these assumptions can be questioned. Indeed, the question of whether or 
not a small wave packet is formed depends crucially on properties of the ground state 
wavefunction, hence on poorly understood features of non-perturbative quantum chro- 
modynamics QCD [|^, ^j. The reduction of the final state interaction is often explained 
as due to the cancellation of gluon emission amplitudes that occurs when a color singlet 
system consists of closely separated quarks and gluons. The legality of adding ampli- 
tudes (before squaring) requires a coherent reaction, so that the cancellation is limited 
to a select class of reactions. An additional issue is the value of the factor M in the time 
dilation; why couldn't it be as large as P? Another way to phrase this question is that 
the expansion of the PLC can be slow only if highly excited intermediate states are not 
important [^, ^. 

The purpose of this paper is to use a simple model with reahstic features to investigate 
the meaning and limitations of each of the three basic assumptions. The model is defined 
and described in Sec. II. Each "nucleon" is the ground state of an electrically neutral 
systems of two quarks interacting via the Coulomb potential. The wave functions for this 
potential can not be computed using pertrubation theory, but can be computed exactly. 
In this sense the above force provides an example of a non-perturbative, yet solvable 
interaction. The two quarks have different masses, with the ratio of the heavier to the 
lighter being two. The nucleus is a collection of such nucleons. This model is relevant 



to color transparency physics because of the analogy between the electric neutrality of 
quantum electrodynamics QED and the color neutrality of QCD. Within this model, it 
is straightforward to construct (Sect. Ill) the wavepacket formed in a high momentum 
transfer reaction. We take the hard interaction as simply adding momentum to the 
heavy quark. This wave packet evolves as it moves through the surrounding nucleus. 
This evolution is governed by the internal Coulomb force between the heavy and light 
quarks and also the final state Coulomb interaction with the surrounding nucleons. This 
interaction is computed, and the dependence on the transverse separation between the 
heavy and light quarks is presented in Sec. IV. The process we consider in this paper 
is meant to be analogous to the (e,e'p) scattering, one of the candidate reactions for 
the observation of CT. In the (e,e'p) reaction hard and soft interactions are of different 
nature - the first one is electromagnetic, the second one is strong. The smallness of 
the soft interaction is a consequence of the color neutrality of the proton, which does 
also have electric charge. In our model the soft interaction is of electromagnetic nature, 
and our "nucleon" is electrically neutral. The high momentum transfer interaction has a 
distinct charge associated with it (only one of the quarks has this charge). The numerical 
results of the calculations are given in Sec. V. The properties of this model are that the 
size of the wave packet created in a high momentum Q transfer process declines with 
Q, and that the expansion rate of this wave packet is inversely proportional to Q. The 
origin of this slow expansion rate is investigated in Sec. VI. A brief summary is presented 
in Sec. VII. 

2 THE MODEL 

We consider a nonrelativistic, electrically neutral system of two quarks interacting with 
the Coulomb potential. The strength of the the electromagnetic interaction, is given by 
the fine structure constant a; a determines the ratio of the binding energy of the system 
to the mass of its constituents. In Nature a = 1/137, but within the present model a can 



be varied and need not be small. Thus the interaction we use is intended to represent 
the strong color force. 

The form of the bound state and continuum wave functions are well known, and are 
not presented here. It is worthwhile to discuss the system of units that we use. Coulomb 
units with h = 1, e = 1, and the mass of the light quark rriq = 1.5 (to set the reduced 
mass rrired to 1) are the natural choice for this problem. With this choice the size of 
the system in the ground state, the Bohr radius tb, is the unit of length, and twice the 
ground state energy is the unit of energy. Note that in the Coulomb units it is the speed 
of light c = l/a that sets the energy scale. In this paper 1/a will be used instead of 
c. A relativistic dispersion relation is used for the energy Ex of the state \X > with 
momentum Px 

El = Ml/a' + Pl/a\ (1) 

The rest mass Mx includes the binding energy ex as well as the quark masses 

Mx = 4.5 + exa'^. (2) 

The energy eigenvalues are 

En = -1/2^2 (3) 

for the discrete part of the spectrum and 

Sk = ^72. (4) 

for the continuous part of the spectrum. 

3 EXPANDING WAVE PACKET 

A wavepacket, possibly a point-like configuration, is created when a photon of momentum 
Q is absorbed by the heavy quark of the "nucleon" in the ground state |1 >. Thus we 
write 

\PLC>=e"^-^^''\l>, (5) 



where spin is ignored, ^f is the position operator for the heavy quark relative to the 
center of mass. Thus Eq. (|^) defines our model of the hard high momentum transfer 
reaction. The transferred momentum is defined as g* = Qz, which specifies the direction 
of the z-axis. 

The use of completeness allows us to express this wavepacket in terms of the elastic 
F{Q'^) and inelastic fx{Q^) form factors 

\PLC >= F{Q')\1 >+J2{fdX)fx{Q')\X > (6) 

X ■' 

and 

F(Q2) = Jd\ri{r)exp{tq- ^rO^i(r), (7) 

fx{Q') =Jd\ri{r)exp{tq- ^f^^x{r). (8) 

where ipx{f) are the Coulomb eigenfunctions (subscript 1 for the ground state). The 
summation (integration) is over the complete set of the Coulomb eigenstates, X = nl for 
the discrete states and X = kl for the continuum states. 

If there were no final state interaction (FSI) of the hadronic wavepacket with the 
nucleus the amplitude of detecting the "nucleon" in the ground state would be equal 
to the form factor F{Q'^). Color transparency is concerned with situations in which 
experimental kinematics constrain the final state interactions to be soft, of low momentum 
transfer. We shall model these soft interactions by treating the nuclear medium here 
as a set of neutral "nucleons". Sect. IV. Then the soft interactions are approximately 
proportional to the product of the wave-packet- "nucleon" forward scattering amplitude, 
with the density of nucleons. In the impulse approximation, the interaction is expressed 
in terms of a matrix element of an operator x{h). The main feature of this operator is 
its dependence on 6, the transverse separation of the quarks (6 ■ 5 = 0). In particular, 
x{h = 0) = 0. We shall discuss a specific model for x i^i Sect. IV below. 

With this notation, the scattering amplitude A^i is given by 

M,{Z)=<l\xG{Z)\PLC>, (9) 



to first order in x where G{Z) is the Green's propagator of a PLC a distance Z (along 
the q direction) through the medium. Color transparency occurs if this term is small 
compared to the Born amplitude F{Q'^). Note that we use lower case letters (b,z) to 
denote transverse and longitudinal quark- ant iquark separations and upper case letters 
(B,Z) to denote the displacement of the center of the wave packet from the center of 
the nucleus. An evaluation of the matrix element of Eq. @ requires an integration over 
cPr, but not over the coordinates (B,Z). (The B-dependence of 7VI1 is suppressed for 
simplicity.) 



In the standard Glauber treatment of final state interactions |Il| an optical potential 
approximation is often used. The potential is proportional to the forward scattering am- 
plitude, hence to the total nucleon-nucleon cross section a. This cross section determines 
the rate of the exponential decay of the scattering nucleon wave function. In the present 
case the b-dependence of \PLC > varies with Z, because the Green's function G includes 
the effects of the heavy quark light quark Hamiltonian. Thus the PLC forward scattering 
amplitude varies with Z. In particular, if the initial state of Eq. (|^) corresponds to one 
of very small transverse size, the PLC expands as Z increases. We shall assume that the 
eikonal approximation is valid. In that case, each of the states from the complete set of 
states defined above in (H), acquires a phase expiPxZ as it propagates a distance Z. 

It is useful to define an effective cross section. aeff{Z), with 

aeff{Z)^M,/F{Q^). (10) 

This quantity depends on the overlap of x(&) with the quark-antiquark wave function. 



and is therefore a measure of how the size of the wave packet varies with Z ||T0[ . Some 
standard manipulations then lead to the result: 

^efAZ)=a+ E i:^#7i^exp.(P„-P)Z + 

/=0,2n=2 ^ \^ ) 

Er*^ilf^-p'(^-^)^- (11) 



This result is similar to the one of Jennings and Miller |T^. In that work the matrix 



elements xx and the inelastic form factors fx were taken from available data and the 
color transparency condition aeff{Z = 0) = was imposed. In the present work we 
use a specific model to evaluate those terms and can determine whether or not the color 
transparency condition is satisfied. 

To proceed, we further specify our notation. The first term of Eq. ([TlD is the total 
cross section for the "nucleon" ground state to interact with a target "nucleon" . This is 

a=<l|x|l>. (12) 

The matrix elements xx are 

Xx=<l\x\X>- (13) 

The orbital angular momentum of the states X are limited to even values by the require- 
ments of parity conservation. We restrict the sum to values of / = 0, 2 to anticipate a 
specific form: x(6) oc 6^. 

The momentum Px of the excited state |X > is given by the energy conservation 
relation imposed by the wave equation: 

P| + Ml/a^ = P^ + Mlla\ P^Q. (14) 

An important quantity is aeff{Z = 0), which measures the size of the initially formed 
PLC. This must be small for aeff{Z) to be small. Note that if aeff{Z = 0) is to be small, 
cancellations in eq. ( pJ]) have to render cre//(0) fall off rapidly with Q ioi Q ^ 1. We 
shall work with a simple form of the soft interaction x{b) = alP' / < 1|6^|1 >; see the next 
section. This form allows the evaluation of creff{Z = 0) with the result 

"'"^'^ = QViTJTT- '^=' 

This is significant, it says that the effective size b^ ~ 1/Q^ for large Q^. This is a 
property of the ground state Coulomb wave function ^ ||. PLC formation is allowed in 
this model. 

For non-zero values of Z the phase factors expi{Px — P)Z spoil the cancellation of 
different terms in the sum {^^. As a result aejf grows as Z increases from 0. This 



indicates that an expanding PLC generally experiences final state interactions. If the 
PLC leaves the nucleus without significant expansion then the final state interactions 
are suppressed. How fast a particular term goes out of phase and upsets the cancelation 
depends on its momentum Px 



with 

j^2 4^2 



AM>^ = ^(4.5-«^^) (17) 



for the discrete spectrum and 

AM^/a^ = (fc2 + 1)(4.5 + aV4(fc2 - 1)) (18) 

for the continuous spectrum. 

Suppose that a limit P ^ oo can be taken, such that all relevant Px -^ P- In 
that case, the Z-dependence of cr^ff disappears and the PLC does not expand. For the 
discrete states P„ ^ AM„/a is true for P as low as ~ S/r^. The situation is different for 
the continuous states, since their energy is not bound from above. The value of P large 
enough to ensure a slow rate of the PLC expansion depends on the structure of the matrix 
elements Xki and form factors fki- This goes back to the idea of Ref. 0, that CT can be 
observed only for momenta transfer greater that the energy of all important intermediate 
excited states. We shall use specific calculations within our model to investigate these 
issues. The use of equations like Eq. (|ll]) in eikonal expressions for PLC wave functions 



is discussed in Ref. |T3|. 



There is another concern about the contribution of the higher excited states. The 
eikonal approximation used to derive the phase factors expi{Px — P)Z breaks down if 
Px -^ 0. Therefore any evidence, that the contribution into the sum (|1T]) of the states 
with Px <^ 1 is important, invalidates the approach developed above. 

We shall also investigate the dependence on a. The physical range of a is from to 
3. The value a = corresponds to the non-relativistic limit in which the speed of light 
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[1/a) is infinite with 



Ex = ex + Tp^- (19) 



The upper hmit a = 3 corresponds to Mi = 0, recall Eq. (^). A further increase in a 
would yield a negative rest-mass of ground state "nucleon" . 

4 THE WAVE PACKET NUCLEON 
INTERACTION 

We are concerned here with deriving the interaction x, which has been defined above 
as the forward scattering amplitude between two quarks (with a mass ratio of two) of 
transverse separation b and a "nucleon" target. We are using a non-relativistic Coulomb 
bound state for the ground state dynamics of the heavy and light quarks, so we also take 
the "nucleon" targets to be the ground state of the same system. 



The expression for the projectile-nuclear forward scattering amplitude is [|ri 

ik 



f{e = 0,b) = — f SB \l - e-^^(^''''^-^)l , (20) 



where the integration is over the area of the nuclear target and k is the momentum of the 
wave packet. The phase shift function xiB-, b,B ■ b) is given by the integral over dZ of 
the sum of the light and heavy quark Coulomb potentials. One usually sees expressions 
in which the integral over b times appropriate wave functions is performed. Here we are 
dealing with a wave packet that is a coherent sum of physical states, so it is convenient 
to study /. Straightforward manipulations lead to the result 

X{B, i. B . j) = Z!£l^ / igi£^ e-f-S (e.i.-.S _ e-J.vS) (21) 

where v is the speed of the wave packet and p{q±) is the Fourier transform of the nucleonic 
charge density, p{r): 

Pir)=^!ilr)-^li^r). (22) 



In particular 



Pili.) = ^—2 - 7 7TTT (23) 

1 + V) 1 + ^ 



where the first (second) term the Fourier transform of the charge density of the hght 
(heavy) quark. The neutrahty condition is 

p{q± = 0) = 0, (24) 

which is vital in obtaining the result that x{B, b,B -b) vanishes at b = 0. This is because 
the integral over cPq± is convergent only because of Eq. (^). Note also that p{q±)/q±_ is 
proportional to the Fourier transform of the Coulomb potential. 

We stress that the expression Eq. ( pT]) is obtained by summing coherently the Coulomb 
interactions of both the heavy and light quarks with target nucleons. This coherence, 
and the consequent "small interactions at small separations" is lost if one is computing 
the cross section for an inclusive process in which one sums over all final states of the 
ejected wave packet. 

We first assess Eq. (^Oj) by expanding the exponent in powers of x- If the usual fine 
structure constant is used, the leading term will dominate. There is no term of O'th 
order in x{B, b,B ■ b). The first order term / d^B x{B, b,B ■ b) vanishes for all values of 
b because the integration over (PB sets g^ to 0. However, / (PB x^{B,b,B ■ b) does not 
vanish. Keeping this term leads to the result 

The zero'th order cylindrical Bessel function has small argument limit 

lim Jo(x) = 1 -xV4 

so that one immediately finds 

f{e = Q,b)r^ib'^. (26) 

This term is purely imaginary, so that its influence is to exponentially damp scattering 
wave functions. 
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We next turn to a complete evaluation of Eq. ( |20D . The integration can be simplified 
by replacing 5 by a shifted value B — ^b so that Eq. (l2ll) becomes 

This expression can be obtained for any ratio of quark masses by using different shifts of 
B. An examination of Eq. ( p7|) shows that x{B, b, B-b) is an odd function oi B-b = cos{(j)). 
The real part of f{6 = 0,b) is proportional to the integral of sin{x{B,b, B ■ b)) over 
between and 27r and therefore vanishes. Thus we have a general theorem that the 
forward scattering amplitude for the scattering of two neutral systems (made of two 
particles) that interact via Coulomb forces is purely imaginary. 

We now evaluate the imaginary part of f{9 = 0, 6). A closed form expression for the 
function x{B, b, B ■ b) can be obtained by performing the integral over d'^q± in Eq. (|27|) . 
The result is 

xiB, b,B-b) = 2a/v[Ko{3x2) - Kq{Qx2) + A^o(6xi) - iro(3xi) + 

1 

-(3x2-ft'i(3x2) — 6x2-ft'i(6x2) + QxiKi{Qxi) — 3xi-ft'i(3xi))], (28) 

where xi =\ B + ^b \ /r^, X2 =\ B — }^b \ /tb and Ki are modified Bessel functions. This 
expression can be used to evaluate the imaginary part of /. The results are shown in the 
Fig. 1. 

We see that one finds a 6^ behavior for small values of b, and in that regime the 
amplitude is proportional to a^. This indicates that perturbation theory is valid at small 
b, even though the coupling constant a can be large. 

Our purpose in this paper is to focus on the properties of wave packets of small 
transverse size. In particular, Eq. (|1^) shows that for large enough Q, the effective value 
of 6^ ~ l/<5^, which is small. Thus we use a simplified version of the interaction 

X{b) ^ fie = 0) ^ ab'/ < b' >, (29) 

where < b^ > represents the ground state expectation value of the operator 6^. 
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The calculations of the present section are a Coulomb version of the calculations 
of Refs. [^, [l^, |18| which used one gluon exchange. Those references also found a 6^ 
behavior for projectiles of small transverse size. The origin is color neutrality, which we 
have modeled here as electrical neutrality. We have made non-perturbative calculations 
of all orders in a which we have taken as large as three. 

5 CALCULATION 

The expression (|11]) for aeff{Z) is evaluated numerically. The matrix elements b"^ and 
form factors fxiQ"^) are calculated for 600 discrete states and a [0,10] range of the 
continuous variable k. This is shown to be sufficient to reproduce the analytic result, 
Eq. (pTSD for momentum transfer Qtb of up to 20. These calculations are performed for 
P = Q which corresponds to the quasielastic kinematics of Ref . |^ . The results for the 
real part of cXeffiZ) for several values of the momentum transfer Q and a = 2 are shown 
in Fig. 2. The "nucleon" expands more slowly (is large for a larger value of Z) for higher 
values of the momentum transfer. There is an initial drop, which is analysed below and 
due primarily to the contribution of the states with / = 2. After the initial drop for 
a wide range of Z the effective cross section ~ Z, which is consistent with the "linear 
expansion" model of Refs. |]T3, |T3[ and also the results of |]T2 . 



The Z-dependence for a given value of Qrs = 6 and three values of a = 1/137, 2, 3, 
is presented in Fig. 3. The figure shows that the "nucleon" expansion distance grows 
with a. This is because Px moves closer to P as a increases, recall eqs. (|TB]). This is 
discussed below in more detail. 

To quantitatively describe the expansion let's introduce an expansion distance Zg^p 
defined as follows 

(^eff{Zexp) = 0". (30) 

The expansion distance Zexp is shown as a function of the momentum transfer Q in Fig.4. 
We see that Zexp is linear with Q; the value of a determines the slope. 
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We calculate aeff{Z) for a reduced range of the continuous variable k G [0, Zcmao;] 
to further investigate the importance of the higher excited states. If kmax is such that 
P{kmax) = kmax ~ VQ then the reduced aeff{Z) is almost identical to the full one. This 
result is important to support the validity of the eikonal approximation. The (Jeff{Z) 
for kmax = 1 is shown in Fig. 5. Even though on the initial stage the expansion picture 
is different, for the most of the range it reproduces the expansion of the "full" PLC well 
and the expansion distance is nearly the same as for the "full" PLC. 

6 DISCUSSION 

A conclusion can be made that a transparency phenomenon is obtained for the model 
under consideration. A small size object is initially formed in a high momentum transfer 
process and the expansion rate of this object is inversely proportional to the momentum 
transfer Q. While the first result has been known for this model for some time [^ Q, 
the second result is a new consequence for this model. We therefore take a closer look 
at what causes such a decrease (favorable for color transparency) in the expansion rate. 
As was mentioned above, PLC expansion depends on the momenta of the intermediate 
states Px- For AMx/ct -C Q ( [iBD can be expanded as 



..^«.^.... (3., 



Thus the values of AMx, recall Eq. (|T^), that correspond to the intermediate states X 
which make important contributions to the sums in Eq. (pA]) determine the expansion 
rate. These states have a discrete or continuum nature. But the energies of the exited 
discrete states are bounded from above and their contributions to aeff{Z); h\i and fni{Q'^) 
decay rapidly and monotonically with the number n. Thus it is more relevant to examine 
the contributions of the continuum states, which in principle can have very high energies. 
The form factors for the continuum fki{Q'^) display a peaking behavior, which occurs 
when the momentum transfer Q matches the relative momentum denoted by the quantum 
number k. See Fig. 6 which shows also that the energy of the states produced in a high 
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momentum transfer process grows linearly with Q. If these states had appreciable matrix 
elements h\i for the soft interaction, the expansion would be very rapid. However, h\i 
decays rapidly with k (Fig. 7). The 6^2 ^^'^ fk2 matrix elements exhibit similar behaviour. 
As a result h\ifki{Q'^) have maxima that experience only slight increase with Q (Fig. 8). 
So we see that in this model high excited states are formed for high momentum transfer, 
but the soft interaction cuts them off. 

The fc— dependence of h^J^i determines the expansion of the PLC. We have seen in 
Figs. 3 and 5 that there is an initial drop for small values of Z . We argue that this drop 
is caused by the contribution of the states with 1 = 2. To see this we expand Eq. (|lT]) 
for small Z: 

XkifkiiQ"^) 






Z' E f dk ^'f^^^^\ e{P, - Pf (32) 

Here we omit, for simplicity, the contribution of the discrete states; kthr is defined by 
^hhr ~ 0- "^^^ calculations show that bl^fko < 0, whereas h\2fk2 > 0. In the first 
integral in ( |32|) the states with / = dominate, which results in the small cre//(0). In 
the second and third integral the situation is reversed, although with the similar result: 
the dominating / = 2 states cause the initial shrinkage of the PLC. 

Another piece of information is provided by the a dependence of the expansion dis- 
tance. The latter decreases with the AMx of the important states \X >. For all the 
discrete states and for the continuum states with k < 1 AMx decrease with a. Since 
Zexp grows with a, it can be seen that only the states with k < 1 are important for the 
PLC expansion. 

This analysis leads to the conclusion that the maximum energy of the states which 
are relevant for the PLC expansion grows much slower than the momentum transfer Q. 
This result is important in two ways. First, it illustrates that this indeed is the condition 
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necessary for CT to occur. Second, our approach was based on the eikonal approximation, 
which assumed the momenta of the propagating states to be much greater than the size 
of the system. This approximation breaks down for the states with energy close to the 
momentum transfer Q. But since these states are not relevant and their contribution 
is negligibly small, the validity of the eikonal approximation is proven by the above 
conclusion. 

7 CONCLUSION 

A "nucleon" model of a hadron has been investigated. In this model hadron consists 
of two quarks bound by the Coulomb potential with the variable strength. A small 
transverse size object is formed when "nucleon" absorbs a high momentum photon. Such 
a system expands with the rate inversely proportional to the momentum transfer Q. This 
slow expansion is a consequence of the fact that the states with energy much greater than 
the ground state energy are not very important. 
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Figure Captions 

• Figure 1. Imaginary part of tlie forward scattering amplitude / for three different 
values oi a/v. 

• Figure 2. PLC expansion for 4 different momenta Q. a = 2. 

• Figure 3. PLC expansion for three values of a. Q = 6. 

• Figure 4. Expansion distance Zg^^p for three values of a. 

• Figure 5. ffe//(^) for the "full" PLC and a "reduced" PLC. Full PLC is constructed 
out of 600 discrete states and [0,10] range of the continuous spectrum. Reduced 
PLC is constructed out of 600 discrete states and [0,1] range of the continuous 
spectrum, a = 2, Q = 6. 

• Figure 6. Inelastic form factors fkoiQ"^) for ten values of Q. 

• Figure 7. 6^q decays rapidly with k. 

• Figure 8. bl^fko for the same ten values of Q as in Fig. 6. 
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